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1. Introduction 

Detailed observations of the properties of the Cosmic Microwave Background (CMB) 
[1] have lifted cosmology to a precision science, raising the bar for theorists who 
compute the predictions with which these observations must be compared in order 
to extract their meaning. This is true in particular for calculations of primordial 
fluctuations from very early inflationary environments, spurring detailed studies of 
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various kinds of corrections to the early classic computations [2]. Such studies [3] — 
[17] are forcing a re-examination of how to interpret the relatively unusual properties 
of fluctuations over super-Hubble distances that have long been known to arise when 
computing for de Sitter, and near-de Sitter, spacetimes [18]. 

Super-Hubble fluctuations raise two, related, puzzles for perturbation theory in 
de Sitter space. One of these is the presence of infrared divergences, and the other 
is the appearance of secular time dependence in successive orders of perturbation 
theory. In this paper we focus on the secular growth, which is troublesome because 
it undermines the validity of perturbative calculations at late times, and so leaves 
open the ultimate late-time fate of the system. We argue that this late-time behav- 
ior can be understood in a controlled way by using the Dynamical Renormalization 
Group (DRG) [19] — a technique for dealing with similar problems in condensed 
matter physics. (See ref. [20] for some earlier applications of DRG methods to cos- 
mology.) The DRG is designed for the situation where perturbative solutions to 
time-dependent equations are limited by the appearance of a secular and growing 
time-dependence that restricts the domain of validity of perturbative methods. The 
DRG allows the late-time behavior to be inferred nonetheless, by using renormal- 
ization group (RG) methods to extend the domain of validity of the perturbative 
solution. This technique is well-adapted to perturbative calculations in de Sitter 
spacetimes, which often reveal such a secular dependence on the cosmological scale- 
factor, a{t). 

We believe it is important to distinguish RG methods, such as these, from whole- 
sale resummation techniques that identify and resum specific infinite classes of graphs 
(such as chain, planar, daisy or cactus graphs). Although RG methods can also be 
interpreted as resummations, these need not be explicitly performed graphically since 
the leading logarithms can instead be obtained by integrating the appropriate RG 
equation. Explicit resummations, such as arise for hard thermal loops or in large- 
expansions for example, are normally required to describe a more serious breakdown 
of perturbation theory than is necessary simply to track large logarithms. 

Although our technique is aimed at understanding the late-time behavior of 
de Sitter correlation functions, we find it also sheds light on the related infrared 
divergences in de Sitter space. This is because the resummed late-time correlations 
can also be less singular in the infrared, and this improved IR behavior can cure the 
infrared singularity. To further explore this connection, we examine the special case of 
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N scalars and compare the DRG resummed results with those obtained using large- iV 
techniques. This comparison shows how the infrared behavior identified by the DRG 
properly captures the effects of dynamical mass generation in the large-iV theory, 
which is also responsible for removing the IR divergences. The DRG resummed 
result resembles a mass even when = 1, however, and this can be interpreted as 
evidence in favour of stochastic arguments for dynamical mass generation for the 
long- wavelength de Sitter modes [21]. However, when we perform the DRG analysis 
for a cubic scalar interaction, we instead find a substantially different result that does 
not have a similar interpretation in terms of a mass. We argue that this is consistent 
with the expectation that a dynamical mass is not generated in this particular case. 

It is important to distinguish between two different types of possible secular 
growth that can appear in de Sitter calculations. The first kind counts the number 
of e-folds between horizon crossing for a given mode with comoving momenta k and 
some conformal time r, generically of the form ln(— A;r). The second kind counts 
the number of e-folds between some fixed time tq, such as the beginning of infiation, 
and the time of interest [22]. For interacting theories the first of these can arise at 
next-to-leading order even in de Sitter invariant situations. The second, by contrast, 
arises in situations when de Sitter invariance is broken by the existence of a preferred 
time scale, and can appear even at leading order without self-interactions. While in 
principle, the DRG could be applied to both cases, in this paper we concentrate only 
on the de Sitter invariant case, tracking loop-induced secular growth from horizon 
crossing. 

Large logs of the form ln(— A;r) have been dealt with in various ways in the lit- 
erature, mostly in the context of observable curvature fiuctuations after a period 
of infiation. In that case, one possibility is to simply keep them in the calculation 
of the scalar field fiuctuations. Indeed for the modes of observational interest this 
factor is at most of order 60 which, while large, may not be so large as to destabilize 
the perturbation series for practical purposes. More importantly, for a finite period 
of infiation the real variable of interest is the curvature perturbation, which is con- 
served outside the horizon for single-field scenarios. In that case one may use the 
quantum calculation to evaluate correlations of scalar field fiuctuations near horizon 
crossing (where the log piece is negligible), then change gauge and follow the classical 
evolution of the curvature perturbation between horizon exit and re-entry [23]. In 
addition, van der Meulen and Smit [5] have argued more formally that the classical 
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equation of motion captures the leading log behavior, up to 1-loop. In this work, 
we use the presence of the ln(— fcr) terms as a tool for understanding the late time 
behavior of fields in de Sitter space. We believe this approach clarifies the physical 
picture in a way that can ultimately be useful for observational questions as well. 

We organize our presentation in the following way. First, §2 describes a repre- 
sentative one-loop calculation for a real scalar field in de Sitter space, illustrating 
the infrared divergences that emerge. This is followed, in §3, by a quick review of 
both traditional and dynamical RG methods, together with their implication for the 
scalar field example of §2. This calculation shows in particular how the late-time 
behavior also resums the large IR logs. §4 then examines the case of scalar fields, 
where it is shown that the DRG reproduces previous results for the late-time limit of 
large- scalar fluctuations. §5 then performs a similar calculation for scalars having 
cubic interactions. We discuss our conclusions in §6. The appendices contain a brief 
review of how the various couplings we examine are renormalized. 

2. Next-to-Leading Order de Sitter calculations 

In this section we use a simple example to display the large logarithms whose inter- 
pretation is the focus of this paper. 

2.1 Scalars in de Sitter space 

To this end consider a single real scalar field, 0, coupled to gravity through the action 



S = — d^x \^—g 



(2.1) 



2 " ■ - ■ 2 " 4! 2 

The subscript '0' is meant to remind that what appears here are bare quantities, 
ml = m\ + 5im? , that include the counter-terms required for the renormalization 
of ultraviolet (UV) divergences. From this point on, we often drop the subscripts 
'i?' and '0,' except when emphasizing the renormalization, and usually quantities 
without a subscript can be assumed to be renormalized. 

Our interest in what follows is in fluctuations about de Sitter space, in which 
case m and ^ only appear through the effective-mass combination^ 

M"^ = m^ + iR = m'^ -UiH"^ . (2.2) 



""^Our metric is mostly plus and we use Weinberg's metric conventions [24], which differ from 
those of MTW [25] only by an overall sign in the definition of the Riemann tensor. Consequently 
for 4D de Sitter i?p„iq„ = H^{gpngmq ~ gpq9rnn), Rmn = SH'^gmn and R = -12H^. With these 
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The massless limit described in subsequent sections describes the case |M^| <^ H^, 
which includes (but is not restricted to) the case ^ = — 0. Notice that this region 
is not an IR attractor of the renormalization group equations for the constants m? 
and ^, which (see Appendix A) over scales H renormalize according to 



(2.3) 



The solutions to these equations are attracted to the case of a conformal scalar — 
= ^ + I = 0, or — > 2H'^ — as // shrinks (provided /i remains larger than H). 
When studying the massless case we must imagine these couplings to lie on an RG 
trajectory that does not get to the attractor before fj, falls below H. 

Scalar fluctuations within this theory may be computed using the 'in-in' formal- 
ism, which aims to evaluate the matrix clement of an operator, 



m 



Texp 



dt'n{t') 



0{t) 



Texp 



dt'n(t') 



m 



(2.4) 

between two 'in' vacuum states. To write this in the path integral form, we double 
the number of fields, {0^,0^}, with one describing the evolution from to t, 
and the other giving the evolution back to tin from t. 
The corresponding correlation functions are 



G-+ix,y)^mx)<f>iy)), (x , y) ^ ^y) <f>ix)) , 

G++{x, y) = e{x' - y')G-+{x, y) + ^(y° - x°)G+- {x, y) , 



(2.5) 



G--{x,y) 
The identity G++ + G~ 



e{x' - y')G+-{x,y) + e{y'- x'^)G-+{x,y) . 

— G^" + G~^ ensures that only three of these are inde- 
K^"*" + 0~) and (pA — (p'^ — (f)~ , this constraint becomes 



pendent. Defining 0^ 
(0a0a) = 0, and in the {(pc, 0a} basis the correlation functions are conveniently 
written 

y \G-+ G- 

where 



(2.6) 



W 



(2.7) 



conventions a conformal scalar has ^ = —h- When necessary, we write the de Sitter geometry using 



spatially flat coordinates and conformal time, ds^ 
and H constant. 



a2(r)(-dr2 + dxidx') with a(r) = -1/{tH) 
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The advanced and retarded propagators are related by GA{x,y) = Gif{y,x) and 
vanish in the coincidence hmit. The free propagators are determined by solving for 
the modes Uk{T) of the field in this background and using 

<P{r, ^) = / (e^'-^«k«k(r) + e-''-aj,<(T)) (2.8) 

where [a^, a]^,] = (27r)^5^(k — k') and «k|in) = 0. Taking |in) to be the Bunch-Davies 
(BD) vacuum, one finds the classic result 



Mr) = Hi'\-kr) (2.9) 

where hI^\z) is a Hankel function, whose order is i/^ = 9/4 — M'^/H'^. For |M^| = 
\m^ — 12^if^| ^ this becomes i/ ~ | — e, with 

-3]p = 3]p-4«. (2.10) 

satisfying |e| <C 1. In the case = — 12^H^ = 0, the modes simplify to 

Mr)^^(l + ikr)e-^'^, (2.11) 



and the Green's functions are 

G°(fc,ri,r2) = -^(G-+ + G+-) 

= ^ I (1 + A;^TiT2) cos[A;(ri - T2)] + k{ri - T2) sm[k{n - T2)] | , 
c^^,{l + 0[{kTr]], (2.12) 
where the last line specializes to the long-wavelength, super-Hubble limit, 

-*r=A«l. (2.13) 

The retarded correlator is similarly 

= e{Ti ~ "^2) -p- |(1 + A;ViT2) sin[A;(Ti - T2)] - k{Ti - T2) cos[k{Ti - T2)] | 

^ 0{r, - T2) ^{tI - r|){l + 0[{kTf]]. (2.14) 

In these expressions the superscript '0' distinguishes the lowest-order, or 'free'. Green's 
function from the loop-corrected one considered in later section. 
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2.2 Divergences 

Some of the divergences arising in loop corrections are already visible when the 
coincidence limit of the two point function is evaluated in real space, 

d^A; 



A(r) = = -iG-+{x,x) = Gl{x,x) 



(27r)= 



Gl{k, T, t) +c.t., (2.15) 



where c.t. denotes the contribution of the counter-terms. The retarded propagator 
does not contribute to this expression since it vanishes at equal times. The isometrics 
of de Sitter space imply (0^(x)) must be independent of x if the expectation value is 
taken within any de Sitter-invariant state (like the BD vacuum in particular) [26]. 

This expression diverges in both the UV and IR when evaluated for a massless 
field in the BD vacuum. To regulate these we introduce both an IR cutoff and a UV 
cutoff, with these cutoffs set (for later convenience) at a fixed physical momentum 
scale: Aj^. < k/a < A^^r (more about this choice below). The regulated expression 
for A(r) then becomes 

a Ain/ 1 7 I / 7 \ 2 



A(r) 



1 



a Air 



^ dk 

~k 



k 



+ c.t. 



(2.16) 
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k 



As before, 'c.t.' denotes the contribution of the counter-terms, an explicit expression 
for which is given in the second line in a particular renormalization scheme, associated 
with a renormalization point fi. We are left with the UV-finite result 

2" 



A(r) 



(2vr)2 
1 

(27r)2 



dk ^ 













Am 



1 + 
1 



k 



(2.17) 



Notice that A(r) is r- independent, as would be expected for a de Sitter invariant 
vacuum, when the divergence is expressed in terms of a physical cutoff. 

As always in quantum field theory, the appearance of an infrared divergence tells 
us something important: since physical observables must be infrared finite, we are not 
yet computing something physical. If our calculation diverges we are not describing 
the relevant long-distance physics with sufficient accuracy, and so must have left out 
contributions whose presence would cancel the IR divergence. But this means we are 
missing contributions that must be just as important as those we compute. 
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Precisely what this missing IR physics is depends on the details of the physical 
question being asked. For instance, if the IR divergence arises in the scattering 
of charged particles, then the missing physics could be related to the possibility 
of radiating soft photons, meaning it was a mistake to choose the initial and final 
states to have a definite number of photons (d la Bloch-Nordsieck [27]). Alternatively 
should the divergence arise when calculating an atomic energy level, it is resolved if 
one uses atomic bound states to perform the calculation. Similarly, for divergences 
arising at finite temperature — either at a critical point [28, 29] or in a hot plasma of 
gauge particles [30] — the divergence indicates the breakdown of the loop expansion. 
In this case loops are not simply counted by powers of the coupling constant, and 
so the missing physics that cures the divergence comes from higher-loop graphs that 
are not suppressed relative to lower-loops by a small parameter. 

For the particular case of a logarithmic divergence, the dependence on the missing 
IR physics is comparatively weak, and more general things can be said about the 
final answer. Once the appropriately modified IR part of the calculation is combined 
with the result above, the logarithms of A^^j must cancel, as in 



to leave Aj^ replaced by whatever the physical scale, L, is that characterizes the 
infrared part of the real physical observable. What is important about logarithmic 
divergences is that this cancellation preserves the coefficient B, so that the coefficient 
of the physical large logarithm, In(yuL), in the full result can be computed purely by 
identifying the coefficient of the IR-divergent logarithm within the UV theory. The 
same is typically not also true for power-law divergences. As applied to atomic 
energy levels this argument leads to the famous factor [31] of In(meL) ~ ln(l/a) in 
the Lamb shift. 

For the case of interest, {(j)'^{x)) in de Sitter space, whether the physical scale L 
depends on time or not depends on whether it breaks the de Sitter invariance. For 
instance, L is r-independent if the relevant IR physics is a scalar mass. Alternatively, 
if it is an earlier pre-inflationary phase that cuts off the long-distance de Sitter 
behavior, then L would in general depend on r. 




(yl + C) + 51n(/iL) + -- - , 
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2.3 Loop corrections 

We now compute the loop corrections, basing our discussion on the one-loop correc- 
tions to the propagator G^, focussing on the contribution from super-Hubble scales. 
The relevant part of the Lagrangian for this theory (after field doubling and changing 
to the basis given by Eq.(2.7)) is 



5c, 
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X 
4! 



) +C.t. 



(2.19) 



Given the Feynman rules from this Lagrangian there are two diagrams that contribute 
at next-to- leading order, one of which is shown in Fig. (1). The second diagram is 
obtained from this one by interchanging ti ^ T2. 




+ 



— X 



1 



Figure 1: The one-loop contribution to the Gc propagator. An additional contribution 
comes from the related diagram with ti and T2 exchanged. Solid lines represent Gc, 
dashed-to-solid lines represent G^, and the crossed circle denotes the counterterms. 



Evaluating Fig. 1 gives the following result for the 1-loop correction [3, 5, 6] 

Gl{k,n,T2) = -^ f drV(r')G°(A;,ri,r')G'°(fc,r',r2)A(r') + (ri^r2) (2.20) 

where the factor of | is from the combinatorics of the contractions, and the quantity 
A(r) is the momentum integral given by Eq. (2.17). The contribution coming from 
long-wavelength, super-Hubble modes is obtained by specializing to that part of 
the integration region for which we can use the asymptotic form for —kTi^2 ^ 1- 
Dropping 0{1) terms relative to ln(— /cr) gives 



Finally, evaluating at equal times, ti = T2 = r, gives the following, IR- divergent, 
next-to-leading contribution to Gc'- 



2 



In ln(-A;r) + --- 



(2.22) 



3(27r)2 

coming from super- Hubble modes. Here the ellipses denote terms whose contributions 
do not diverge as either Aj^^ or —kr vanish, and so are subdominant to those explicitly 
displayed. 

2.4 Large logs and secular time dependence 

Eq. (2.22) displays two separate kinds of large logarithms. The one of interest for 
the arguments of the next section is the ln(— /cr) term, since this displays the secular 
time-dependence for which the dynamical renormalization group is designed. This 
term is proportional to the number of e-folds between horizon exit of the mode k and 
the time r where we evaluate the result: —kr = k/{aH) = ak/a. For calculations of 
curvature perturbations resulting from inflation, this apparent growth does not give 
rise to a physical divergence: the curvature perturbation is conserved outside the 
horizon and in fact one may follow the classical evolution of the modes shortly after 
horizon crossing until horizon re-entry [23]. 

The other large log is the IR divergence as Aj^ 0. This has its origin both 
in the use of massless fields, and in the eternity of de Sitter space that the above 
calculation implicitly assumes holds in the past of r. As discussed earlier, this 
divergence ultimately cancels some other dependence on in a physical observable, 
and because the divergence is logarithmic the cancellation of A^^ gives 



2 



A 



In (/iL) In {-kr) 



(2.23) 



3(27r)2 

where L is a physical scale associated with whatever long-distance physics makes 
the result IR finite. The precise nature of this scale is model-dependent, because it 
cannot be identified purely within the UV part of the theory we use to this point. 
We identify below what L is in the specific case where it is a mass that provides the 
convergence. 

3. RG Methods 

This section briefly recaps standard lowest-order RG arguments, for their later use 
in our cosmological applications. Part of our purpose is also to emphasize how the 
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use of the RG equations avoids the necessity of performing exphcit summations over 
infinite classes of graphs. 

3.1 The standard argument 

Standard one-loop calculations for dimensionless couplings give the value of the renor- 
malized couplings, for different values of the renormalization point, fi: 

a(fi) = a(fio) + 6a^(/io) In ( — ) , (3.1) 

where 6 is a calculable number. The domain of validity of this calculation requires 
both a{fio) ^ 1 and a(/io) ln(/i//io) <^ 1. 

The RG proceeds [29] by differentiating the above result with respect to fi, giving 

figj^ = ba\ (3.2) 

where, to within the accuracy given, a on both sides can be taken to be a(/i). This 
equation can be integrated to give the solution 

— = — -bln(-^ 
a(/i) a(/io) V/^o 

The point of this seemingly content-free exercise of differentiating and then integrat- 
ing the loop result is this: while eq. (3.1) relies on the product a ln(/i//io) being small, 
eqs. (3.2) and (3.3) only require a be small and so have a larger domain of validity, 
applying in particular when aln(/i//io) is order unity. Integrating the RG equation 
effectively resums the leading logs. Notice that this resummation selectively takes 
terms to all orders in a ln(/i/ /io), but this does not rely on identifying which subclass 
of graphs is dominant. 

3.2 The dynamical RG 

The dynamical RG [19] uses a similar logic to identify the long-time behavior of 
solutions when perturbative methods generate perturbations that grow secularly with 
time. For instance, suppose an approximation scheme based on expanding in a small 
parameter e produces a result of the form 




y(t) =yo(t)+eyi{t) + 0{e^) 

= yo{c,t)+eyi{c,t) + 0{e^) , 



(3.4) 



where the second Une emphasizes the dependence of the solution at each order in 
e on the integration constant, c, of yo{t). Of particular interest is when yi{t) is a 
secular, growing function of time. When this is so, the expansion in powers of e 
breaks down for time-scales, t>T, where eyi{T)/yo{T) is order unity. 

The dynamical RG shows how to infer the behavior of the solution for times 
t > T. To do so, first introduce an arbitrary time scale, with 



y{t) = yo{t) + e [yi{t) - yi{d) + y,{d)\ + O {e^) , 
and absorb the last yi{'&) term into the zeroeth-order integration constants: 

yo[c(i?),t] :=yo(c,t) + £yi(^?), 

and so 



(3.5) 



(3.6) 



(3.7) 



y{t) = yo[c(i?), t] + e [yi{t) - y^{^)\ + O {e') . 

The dynamical RG argument proceeds from the recognition that y{t) is indepen- 
dent of dy/d'& — 0. Differentiating and dropping 0{e'^) terms then gives 



dyo 
dc 



dc _ dyi{c, ^) ^ 



(3.8) 



which can be regarded as a differential equation to be solved for c('t?), whose solution 
is c = c{{}). 

The argument is now the usual one: the solution c = £(■(?) has a broader domain 
of validity (typically simply e <^ 1) than does the initial definition of c(^?) (which 
required eyi{t) <^ 1). The final late-time behavior is then found by using the fact 
that y{t) does not depend on d to make the convenient choice d — t, and so 



y{t) = yo[m,t] + e [y^{t) - y,{^)\ + O {e') 
^yom,t] + 0{e') . 



(3.9) 



This result resums the leading large-time behavior beyond the leading order pertur- 
bative result. Particularly simple is the case where 



y{t)^c [l + sf(t) + 0{e')] , 



(3.10) 



where yo = c is time-independent. In this case 0(1?) = c[l + e /{-&)] and so the 
condition dy/d-d = implies 

r\r H f 

(3.11) 



dc 
d^ 
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which integrates to give c = c('i9), where c{'&) = ce'^^^'^\ This gives the resummed 
large-time solution 

= ce^^W [l + 0{e^)] . (3.12) 

This is the particular case that arises in the one-loop scalar calculation on de Sitter 
space. 

3.3 The dynamical RG for de Sitter fluctuations 

Earlier sections show the leading infrared-sensitive part of the de Sitter space calcu- 
lation of massless-scalar fluctuations involves a secular time dependence, 



2A;3 



1 



A 



In (/iL) In {-kr) H 



(3.13) 



3(27r)2 

with ellipses indicating terms that do not diverge when L or —kr vanish. For this 
result, and assuming that the scale L is independent of t — as is the case if the 
IR physics is de Sitter invariant (like a scalar mass term) — the dynamical RG 
resummation described above predicts the following long-time behavior 



Gc{k,T) 



2k^ 

El 



exp 



with 



k 



5 



A 

3(27r)2 

5 



In (fiL) In {—kr) 



1 + 



1 + 0{6^ 



A 



In (fiL) 



(3.14) 



(3.15) 



3(27r)2 

This resummed result holds to all orders in (51n(— A;r), but neglects terms that are 
suppressed by additional powers of 6 without accompanying powers of ln(— fcr). Its 
validity therefore requires 5 -C 1, although this condition might also be relaxed 
through an RG resummation of the large In (fJ^L) logarithm. 

Notice that because fiL ^ 1 and A > the power 6 is also positive, implying 
that the DRG-improved correlator, eq. (3.14), is less singular as — than is the 
zeroth-order result, G^{k,T) oc H^/k^. Also notice that although we do not need to 
know which subclass of graphs is dominant in making this argument, the results of 
ref. [6] show that these leading logs come from the 'chain' diagrams that repeatedly 
insert a scalar self-energy. The validity of the DRG does not require that these 
entire diagrams dominate all of the others, but only that they dominate in their 
contribution to the leading logs. 
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4. Comparisons to other results 

This section makes three points. First, the small-A; hmit of the above DRG expression 
for Gc{k, t) is argued to be similar to the small- A; result for a massive scalar, showing 
that the resummed late-time behavior softens the IR divergence in way very similar 
to a mass. Second, the super-Hubble part of the loop calculation is repeated using 
massive scalars, to show precisely what the large logarithm ln(/xL) means in an 
explicit example for which the IR physics is known and calculable. Finally, we 
extend our calculations to an A^-scalar model for which a dynamical mass is known 
to be generated within a controllable approximation in the large- limit. We do so 
in order to show more precisely how the DRG resummation reproduces the effects of 
the dynamically resummed mass. 

4.1 Massive scalar field 

A massive scalar field in de Sitter space provides the simplest example of IR- convergent, 
de Sitter invariant, physics which nonetheless can have interesting contributions from 
super-Hubble modes if the mass is less than the Hubble scale. Just as is true for flat 
space, a mass term improves the IR convergence by making propagators less singular 
for small k. To see why this is so for de Sitter space, recall that the massive scalar 
propagator obtained from the general mode functions, eq.(2.9), becomes, in the limit 
-kr = k/{aH) < 1 

G0(fc,n,r2) ^ ^s(^'r,T,Y (4.1) 

G°(A;,n,r2) ^ e{n - r,)^{rt'T^ - r^-) . (4.2) 

where e = M'^/3H^. As advertised, these expressions cure the IR divergences en- 
countered previously because they are less singular than the massless case as A; ^ 0. 
But for <^ the growth of the difference between the massive and massless 
expressions for Gc{k, r) is much slower than in flat space, with significant deviations 
only arising once {—kr)'^'^ deviates from unity. This occurs for k < k^, with 

aH 

and so k^:/a ~ }j ^-^^'^/'^^i'^ can be much smaller than M. 

What is important for the present purposes is that the small-/c form of the mas- 
sive propagator has the same kind of /c-dependence as does the DRG-resummed result 
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for the massless field: both behave as Gq oc /c""'^, for a > 0. For massive particles 
we have a = 2e = 2M'^/3H^ while DRG resummation gives a = A ln(yuL)/127r^. 
Equating these powers shows that the DRG-resummed long-wavelength correlator 
behaves as if it describes a dynamically generated mass of order 

M^^ ^ ^ ln(/iL) . (4.4) 

A reasonable guess for the size to take for L in this expression is 

L ^ ^ ^ i eV2^ , (4.5) 

(a more precise choice for L is described shortly), in which case we have 

4.2 Massive loop calculation 

It is instructive to repeat the IR-sensitive, super-Hubble part of the loop calculation 
using this same massive scalar (still with a quartic self-interaction), since this provides 
a simple and specific example of IR-safe physics that can cure the IR divergences 
of the massless case. Among other things, such a calculation allows us to identify 
more precisely what the large logarithm ln(/iL) means when it is a small mass that 
is responsible for IR convergence. 

So far as the purpose of identifying the IR behavior is concerned, what is impor- 
tant is following the contribution of the super-Hubble modes in the result. For the 
loop factor, A(r), the contribution of these modes is 

(27r)2 J^. p (27r)2 J^. p 



2e V27r / 



IR ^ \ ) -J a i^m 

/iX2e /A,«,^'^^ 



HJ \ H 



(4.7) 



which is again r-independent because of the choice that and /x are physical 
scales. Notice that in the limit e — eq. (4.7) approaches the super-Hubble part of 
the massless case, eq. (2.17), 



ZT \ 2 / 



but instead converges in the IR if e is held fixed: 

lim A = — I — I , 

Affi^o 2e\2'n) \H 



- 15 - 



Here /x is, as before, a renormalization scale associated with the UV counterterms. 
Using eq. (4.9) to evaluate the 1-loop result. Fig. 1, then gives 

GUK ri, r2) = -\j dr'a\T') r') G° (A:, r', r^) A{t') + (n ^ r^) 



A 
2 



1 /i/ 



2e \27r 



2e 



dr' 



6A;3 



-i/fc 



3-e 



r' 

in ^ T2) 



4e V27r J \^ 



(4.10) 



where the unwritten terms are finite as /cr ^ 0. We find in this way the loop- 
corrected result 



Gc(fc,r)^^(-A:r) 



2e 



A 



6(27r)2e 



(|)'^n(-A;r) + --- 



(4.11) 



Notice that this is similar to expression (3.13) in the massless case, with the replace- 
ment 

For small e this reduces to 



(4.12) 



ln(/iL) ^ ^ + In (I ) + 0(e^ 



(4.13) 



which agrees with the simpler estimate, eq. (4.5), of the previous section. In this case 
the large logarithm In(yuL) is seen to correspond in the full result to the enhancement 
of the loop contribution by the factor 1/e. 



DRG resummation 

Applying the DRG to resum the leading logs gives in this case the late-time result 



Gc{k,T) 



2A;3 



-kr) 



2e+5„ 



(4.14) 



where 



A / \ 



2e 



confirming that the resummed leading logs modify Gc{k, r) in the same way as would 
a mass. Notice that although this resums terms to all order in 6^ ln(— A;r), it requires 
5m 1, and so breaks down once A/e becomes too small. When both masses and 
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couplings are present the mass effect dominates the couphng effect so long as 5^ < 2e, 
or 

For A larger than this it is the coupling that dominates in the IR, and it contributes 
in the same way as would a mass of order 

,2 SH^ , A/72 /^N2e 3^^4 



where the last approximation neglects corrections that are of order e, since the as- 
sumption that the coupling effect dominates requires to be systematically small 
relative to A/ (47r)2. 

Notice that eq. (4.17) precisely agrees with what would be expected in a mean- 
field approximation, for which^ ^mf — k -^('^^)' given the nonzero expectation 

that a massive scalar field prepared in the Bunch-Davies vacuum acquires in de 
Sitter space [34, 35]. To this extent the DRG-improved long-time limit agrees with 
what would be expected from a stochastic approach [21] to super-Hubble scalar 
fluctuations. We explore this in more detail using a large- limit [21, 10, 6] in the 
next section. 

4.3 Dynamical mass at large 

Large- expansions provide a useful laboratory for understanding dynamical resum- 
mations within a controlled approximation, and so we therefore pause to examine 
them briefly here. In this section we apply the dynamical RG to the large- general- 
ization of the above calculation in order to show how it compares with the long-time 
behavior predicted by the DRG. 

Consider to this end the large- A^ generalization of our model. 



(4.19) 

^For the numerical factor, compare with the large- iV expression in the next section. 
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where = — 12^if^ in de Sitter space, $ denotes an A^-component column 
vector of real scalar fields, and the first line follows from the second if s is integrated 
out using the exact result 

s = $'^<l> . (4.20) 
The scalar field equation for $ in this model is 

+ I^M^ + $ = 0. (4.21) 

In the limit of large we take A small and (s) large, so that g = NX and cr = (s) /N 
are fixed. In this large- iV limit the field equation behaves sensibly, taking the mean- 
field form [36] 



+ Af 2 + — $ = , (4.22) 
6 J 

which uses the large- expression As $ = (^cr $ (1 + 2/A^) + ■ ■ ■ . This shows that the 
scalar self-interaction contributes to the mass for $ in the large- iV limit when terms 
of order are dropped 

Ml = + ^ = - 12^H^ + ^ , (4.23) 
6 6 

provided only that a is nonzero. 

The static value for a can also be computed for large in de Sitter space, with 
a standard result given by the contribution of A^ free scalar fields [34, 35]: 

and so 

M^ = M^ + ^ = M^ + ^J^. (4.25) 

Using this mean-field mass in the $c correlator then gives the large- A^ expectation 
for the form of the small-/c limit of the de Sitter correlator 

Ge(A;,r)^^(-A;r)2^-, (4.26) 



where 



2Mi 2M' gH^ 

3H^ 3H^ 6(27r)2M2 ' ^ ^ ^ 
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DRG resummation 

For comparison we instead compute Gcik, r) perturbatively in and repeat the loop 
calculation of the previous section in the large- iV limit. There are only two changes 
required. One change comes from evaluating the combinatorial factors arising from 
performing the contractions of the various fields appearing in Fig. 1, which converts 
a factor of 3 (from counting the number of ways to contract fields) in the previous 
loop calculation into (2 + A^). The second change is the replacement of the coupling, 
A g/N, and so the generalization to arbitrary of the previously obtained loop 
result is obtained by making the simple replacement 

3\^g(l + ^] . (4.28) 



This leads to the following large- generalization of eq. (4.11) 



1+ (^) '\n(-kT) + 

18(27r)2e \hJ ^ ' 



(4.29) 



Using the DRG to resum the large-r behavior gives, as before 

Ge(A;,r)^^(-A;r)2^+^-, (4.30) 

where 

5 =^ (4 31) 

18(27r)2eWy " 6(27r)2M2 ' ^^""'^ 

and the second, approximate, expression uses e ^ 1 to drop the factor [fi/H)'^'^ 
(compare with eq. (4.15)) . For small k the total power of k therefore is 

2M\ gH^ 
2e + 5. = — + ^^^^, (4.32) 

in exact agreement with e,^ obtained using the large- A^, mean-field result, eq. (4.27). 
Self-consistent masses 

The large- A^ limit also allows a more careful exploration of the case of very small M^. 
Because the dynamical mass, of eq. (4.25), is bounded from below, the physical 
mass does not disappear even in the limit — ^ 0. Instead it reaches its minimum 
value when OM^/dM'^ = 0, or 

Ml^ = ^^H\ (4.33) 
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Because g/^ — ^ for this value of it lies within the domain of validity of our 
approximate expressions, which require gje ^ 1. The subsequent growth of for 
smaller than this is suspect because gje becomes larger than unity. This means 
that eq. (4.33) corresponds to the smallest mass that it is possible to get within 
the domain of validity of our approximations. (The corresponding result for = 1 
would be M^jj^ = -\/3A H"^ /An. This mass is one that has been obtained by solving 
a 'gap' equation, in ref. [6], and agrees with the much earlier result of Starobinsky 
and Yokoyama [21], after rescaling to match their definition of A.) 

5. Scalar field with a cubic interaction 

In this section we consider a slightly different scalar field with a purely cubic 

interaction 

V{<i>) = j,<t>\ (5.1) 

with no mass term and no quartic interaction. We find in this case that the DRG- 
improved super-Hubble contributions do not have a momentum dependence appro- 
priate to a dynamically generated mass. 

-o — o — o- 

(a) (b) (c) 

Figure 2: The three contributions to Gc(A;,ri,T2) at one loop. Diagram (a) is divergent 
in the UV and IR, (b) is only IR divergent, and (c) diverges only in the UV. 



The one- loop graphs (other than the counterterms) contributing to Gcik, ti,T2) 
in this theory are shown in Figure 2. We evaluate these graphs^ using only the long- 
wavelength, super-Hubble limit of the propagators, as before. All graphs contribute 
a factor of [ln(— A;r)]^ thanks to the two interaction vertices, but each has a differ- 
ent structure of divergences in the momentum integral. The first diagram gives a 
contribution that grows in time and is proportional to A = (if/27r)^ ln(yu/A^fl), as in 

^The complete Feynman rules for this theory in the real time formalism are given in [5]. The 
left vertex of Figure 2(a) contributes — i/ia'*(ri)(5(ri — t')5{ti — t"), where ri is the time associated 
with the incoming line and t' ,t" are associated with the lines to the right of the vertex. This is 
the only new rule needed to evaluate the IR divergent contributions in this example. 
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the (j)^ case. The second diagram converges in the UV and the momentum integral 
is dominated by loop momenta p ^ k. The last graph is UV divergent but not IR 
singular, and so does not contribute to the terms of present interest. 

Evaluating the contribution of the first and second graphs gives the equal-time 
result 



Ga{k,T,T) = GUk,T,T)\l + 



h 



2 r 



9H^ 



\n\-kT) + -ln(-A;r) 



(5.2) 



+ -{-kT)^\n{-kT) 



4A + 



In 



aA 



IR 



G°(fc,r,r) 1 



9(27r)2/f2 



\n^{-kT) + 



~9lF 



\n\-kT) + ... 



where A(r) = (/J/27r)2 ln(/i/A jr) is as defined in the previous section, and the second 
line uses k/{aAjji) = (— A;r)(/J/A//j). A massive scalar field with a cubic interaction 
is treated in detail in [5], who also do a careful analysis of the UV pieces. However, 
the IR part of the calculation above is sufficient to contrast the result of the DRG 
for this case with the quartic interaction result. 

Applying the DRG method to this expression exponentiates the quantity in 
parenthesis. 



Gc{k,T,r) = G^^{k,T,T)exp 



9H^ 



J^,^A-kr) + f,ln\-kr) 



(5.3) 



which clearly does not share the small-k dependence appropriate to a constant mass. 
This agrees with naive expectations, since even if a mean-field treatment were valid 
it would predict M^^^ oc h{(f)) = 0. But unlike the quartic case, for the purely cubic 



theory the semiclassical expansion about the solution (</>) 
the study of its fluctuations problematic. 



is unstable, making 



6. Conclusions 

Loop calculations involving massless scalars in de Sitter space have long been known 
to be plagued by strong infrared effects. These arise due to the indefinite growth of 
fluctuations on super-Hubble scales. These long-wavelength effects complicate the 
precision calculation of inflationary observables, by potentially putting the late-time 
behavior beyond the reach of perturbative calculations. 
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We argue here that the DRG is a natural tool to use for this kind of problem, 
where secular time dependence limits the domain of validity of a perturbation ex- 
pansion. It provides a simple way to resum these terms, thus allowing the late-time 
behavior to be extracted from perturbative calculations. 

We apply the DRG to several simple examples, such as a single quartically self- 
coupled field as well as a large N quartically coupled model. In these situations 
the resummed IR behavior is less singular than is the leading-order result, in a way 
that mimics the presence of a nonzero mass. The same analysis for a cubic scalar 
theory demonstrates that the cubic interaction apparently does not self-regulate by 
generating a mass. 

Fortified by these checks, we believe that DRG is likely to prove to be very useful 
for more detailed studies of other puzzling aspects of late-time, long wavelength de 
Sitter fluctuations. 

Another direction to take this work might be to consider situations where there 
the important long-distance physics breaks de Sitter invariance. This could happen 
because inflation is not eternal, instead beginning at the end of a pre-infiationary 
phase [22], or perhaps if dS space should prove to be only metastable, with the BD 
vacuum ultimately decaying on large scales due to the inhomogeneous structure, 
perhaps coming from the formation of bubbles of some new phase. Whatever this 
source of dS symmetry breaking is, the two point function becomes a function of 
time and one may wonder whether a DRG analysis can still be useful. 

An interesting application of this type would be to reproduce past discussions of 
such secular behaviour, such as that of ref. [37] which discusses large infrared loga- 
rithms in a time-dependent situation, and argues, following a proposal of Starobinsky 
and Yokoyama [21], that these may be resummed [38] using ideas from stochastic 
inflation. Investigations of this point, and how it relates to the DRG techniques 
discussed here, are in progress. Previous arguments about applying RG techniques 
for the case where the dS symmetry is broken can be found in [40]. 
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A. Renormalization 

When renormalizing the one-loop UV divergences it is important to keep all of the 
interactions whose counter-terms can be required to cancel divergences at one loop. 
For an 0(A^)-invariant system of scalar fields the required terms are: 



S = - d^'x 



2 g'-'^d.^^d,^ + + ^R){<^^<^) + ^ (A.l) 
+ Vo + AR + B R^^xpR"''^^ + C Rf.^R'"' + DR^ + EUR 



The coefficients of these terms contain divergent contributions, whose presence can- 
cels all UV divergences that arise at one loop. The counterterm lagrangian consists 
of the same terms with —>■ 5m?, C, —>■ and so on. 

Working in dimensional regularization, a standard calculation [39] gives the di- 
vergent part of the one-loop effective action, Sqo, generated by an arbitrary collection 
of fields. Assuming these fields have an action of generic form 

S=-\j ^^^V^ $^(-n + X)$ , (A.2) 

one finds^ 



Soo = ±2 (^^J T{2-n/2)fi^-^ J d^xV^tTa2, (A.3) 

where the upper (lower) sign applies for boson (fermions), n —>■ 4 denotes the dimen- 
sion of spacetime, fi is the arbitrary mass scale that arises in dimensional regulariza- 
tion, and a2 is the second Gilkey/de Witt coefficient, 

a2 = ^ (2R^,^pR>^'''' - 2R^,R^'' + 5R' - UDR^ 

+ ^RX+^-X^-^nX + ^Y^.uY^'' . (A.4) 



^This can be written in an alternative form also involving the lower-order Gilkey coefficients ag 
and ai by grouping terms involving different powers of m^. 
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Here [D^,D,^]^ = Y^y(^ defines Y^y, where is the covariant derivative appearing 
within □$ = g^"'Df,Dy^. 

Since gauge- neutral scalars are the case of present interest, we have Y^^ = 0. 
Further speciahzing to the potential V = \{m? + ^ + A($^$)^ then gives 



A 



6 



and so 



and 



A 



tr X = Nim^ + iR) + -($^$) (A^ + 2) 

6 



trX^ = N{m' + iRf + - (m^ + iR){^^^) (iV + 2) + —{^^^f {N + 8) 

3 36 

Using these in eq. (A. 4) gives the result 

tra2 = ^(2i?^,Api?^"^'' - 2R^yR^'' + hR^ - 12UR^ 
+ ^i?trX + ^trX'- ^ DtrX 



(A.5) 



(A.6) 



(A.7) 



— ( E ^ B^^^P - H E^" 



A^ 



A 



6 



R' 



N 
6" 



UR 



N[i + -] m'R + -(A^ + 2) m^($^<l>) 



6 



6 
36 



A^ 



+ ^(A^ + 2)[i + ^]R ($^$) + ^(A^ + 8)($^$)' 



6 



(A^ + 2)n(<i>^$). 



(Ai 



Using the above results to compute Sqo, shows that the one-loop divergences are 
given by poles as n — > 4 of the form 



"/2 / Ti-4 



n - 4 



d'^xv^^ tr a2 . 



(A.9) 



These can be absorbed into minimal-subtraction (MS) counter-terms of the form 



■'ct 



1 

47r 



d'^Xi/^ 



n- A 



In 



tr a2 



(A.IO) 



and the counter-terms in other renormalization schemes differ from these by finite 
renormalizations. Notice that for the counterterm action unitarity requires that all 
fields are evaluated at n = 4. 
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Comparing this with the counter-term action and differentiating with respect 
to II shows (in minimal subtraction) that the renormahzed couphngs satisfy the 
renormahzation group equations 







dji 


487r2 


dX 


X^N 




~ 487r2 


dA 


w?N 






dD 


N 




327r2 



1 + 

1 + 



_ XN 
dVo m^N 
dB 
dE 



1 + 



N 



327r2 
dC 



N 



dn 28 807r2 



N 



967r2 



1 



(A.11) 



These may be integrated in the usual way. Write the first few of these equations 



as 



and 



where 



6 

h X2 



dm? 
djj, 



— biX , 



&i = 



N 
48^ 



1 + 



N 



The A equation integrates to give 

A(/.) = 

and the other two give 



and bo 



An 



N 
487r2 



TV 



1 - 62 Ao ln(/x///o) ' 



(A.12) 
(A. 13) 
(A.14) 

(A.15) 



where 



e(/x) + 1/6 _ m\fi) 
Co + 1/6 ml 

61 A^ + 2 



An 



^ 62 A^ + 8 



= 1 



6 

TV 



+ 



(A.16) 



(A.17) 



and so 1 > p > |, with the smallest value occurring when A^ = 1. 
In particular, 

A\^ 1\/A^^ 



- 12 = ml 



Ao 



12 Uo + 



6 



Ao 



(A.18) 



Recalling that both 62 and p — 61/62 are positive, we know A is not asymptotically 
free, and so // < //q imphes both A < Aq and (A/Aq)^ < 1. Consequently as A flows 
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to and ^ flows to — | with falling /i, the quantity — 12 ^H"^ instead flows to 
2H^ > 0, ensuring ^ is necessarily much smaller than 12 \C,\H^ 2H^. 

The Riemann- and Ricci-squared coefficients are also easily integrated, giving 



B = Bo + 



N 



28 807r2 



In 



and C = Co 



N 



In 



ii 



(A.19) 



28807r2 

Those for A, Vq and D require integrating the previous expressions for m^{ii) and 
{(/x), with 

N 

2p 



A{n) = >lo + 



167r2 



16^^ 



= ^0 + 



167r2 





rA(/x')i 







A^mg 



167r2Ao(26i - 62) 



ln(////io) 



6 + 



dt 



327r2 



K)(/^o) + 



A^m^ 



e; Ul - Mo ln(////xo)]^^'-^ 



- 1 



Ao 



(A.20) 



327r2Ao(26i - 62) I [1 - b2\o Hi^/i^oWp-^ 



N 
32^ 



Co 



N 



6 



^d/i' 



Ao 



2p 



e; I [1 - ^'2Ao ln(/x//Xo)]2f-i 



- 1 



327r2Ao(26i - 62) 
Running in the large- A^ limit 

When taking the large- A^ limit we write X = g/N and may take 61 = 62 = Nb, up to 
sub dominant terms in 1/A^, with 

1 



487r2 



(A.21) 



Because bi = 62 at large A^ we also have p = 61/62 = 1 in this limit. The first few 
RG equations in this limit then become 



bgm^, fi 



bg{^ + 



dpi ^ ' "^dpt 
and so on. These integrate to give 

+ 1/6 m\pi) g{ii) 



9g , 2 



Co + 1/6 



rrin 



go 



1 - 62 5'o ln(/////o) ' 



(A.22) 



(A.23) 
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which show that m^, ^ + | and g all renormalize in the same way, up to cor- 
rections. This implies that the relative size of all of the terms in expressions like 



do not change as /i scales, and and both approach 2H^ as /i ^ 0. 
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